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Abstract 

Recent developments in the theory of heavy quarks are reviewed. In the area of 
heavy quark fragmentation, there has been progress in the study of both pertubative and 
nonperturbative processes, including the identification of new observable nonperturbative 
fragmentation parameters. There has also been considerable theoretical activity in the 
study of inclusive rare and semileptonic decays of bottom hadrons. 
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1. Introduction 

The properties of hadrons containing a heavy quark Q simpHfy considerably in the 
hmit niQ ^ oo. If the hadron contains only heavy quarks, this is because of asymptotic 
freedom. If the hadron contains light quarks, this is because for niQ ^ Aqcd the light 
degrees of freedom become insensitive to the mass mg; as far as they are concerned, the 
heavy quark acts simply as a non- recoiling source of color. In the latter case, hyperfine 
effects associated with the heavy quark chromomagnetic moment also decouple, and a new 
"heavy quark spin-flavor symmetry" emerges. There is an extensive literature in which 
these simplifications have been used to predict hadronic spectra and decay rates, among 
other things. In this talk I will describe more recent work, in which they have been applied 
to heavy quark fragmentation and inclusive decays. 

2. Heavy Quark Fragmentation 

There has been much recent progress in describing the fragmentation to hadrons of 
heavy quarks which are produced at high energies in collisions. In the case of fragmentation 
to a hadron containing a single heavy quark, the symmetries of the heavy quark effective 
theory have been applied to identify new nonperturbative fragmentation parameters which 
describe the anisotropies of the light degrees of freedom. For fragmentation to hadrons 
containing two heavy quarks, such as to "onium" states, perturbative QCD has been used 
to compute fragmentation functions and probabilities in an expansion in asimq) and 
mq/E, where E is the energy of the quarkonium in the rest frame of the event. We will 
discuss each of these systems in turn. 

2.1. Heavy-Light Systems 

We begin by applying the heavy quark symmetries to fragmentation processes in which 
no new heavy quarks are produced.^^ In the limit mg oo such a process factorizes into 
short-distance and long-distance pieces. The heavy quark is first produced via some high 
energy interaction, such as the decay of a virtual photon, which may typically be calculated 
in perturbation theory. This initial stage occurs on a time short compared to the time 
scale 1/Aqcd of the nonperturbative strong interactions. Later, over a longer, hadronic 
time scale, a fragmentation process occurs which eventually produces a physical hadron 
containing the heavy quark. Since this slower process only involves the redistribution of 
energies small compared to mg, the velocity of the heavy quark remains unchanged once it 
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has been produced perturbatively, and its mass and spin decouple from the nonperturbative 
dynamics. 

Let us now consider this sequence of events in somewhat more detail. We imagine 
that we begin with a heavy quark which has been ejected at relativistic speed from a hard 
reaction. The axis linking the rest frame of the heavy quark to the center of mass frame 
of the hard process is a preferred direction, which we call the axis of fragmentation. The 
interactions which couple to the heavy quark spin, such as the chromomagnetic moment 
operator, are suppressed by at least one power of the small parameter Aqcd/wq, and 
hence the rate of heavy quark spin flip is very slow on the hadronic time scale 1/Aqcd- 
We might imagine the early stages of fragmentation to involve the production of excited 
hadrons containing the heavy quark, which then rapidly decay to lighter excited states. 
So long as these decay times are shorter than the typical spin flip times 'tt-q/Aqqj^, the 
heavy quark spin will be unaffected by this evolution. 

However, eventually the light degrees of freedom will reach a state whose lifetime is 
longer than the time required to flip the heavy quark spin. In the heavy quark limit, if 
the light degrees of freedom have angular momentum j > 0, the physical hadrons are a 
nearly degenerate doublet with spin J = j ± ^. The trademark of a sufficiently long-lived 
state is that the members of this doublet are well separated resonances; in that case, the 
doublet lives long enough for the two eigenstates of total angular momentum to become 
incoherent with each other. This, in turn, will happen only if the heavy quark and the 
light degrees of freedom have enough time to exchange their spin orientations before the 
doublet decays. 

It is at this point that things become interesting, since the dynamics now allows the 
different helicity states of the light degrees of freedom (and of the heavy quark) to mix with 
each other. This effect, as we shall see below, is observable in the subsequent decays of the 
long-lived states, in such a way so as to give us information about the fragmentation process 
itself. We now introduce a new set of fragmentation parameters appropriate to HQET. 
When light degrees of freedom with angular momentum j are formed by fragmentation, 
they can be produced in one of 2j + 1 helicity states along the fragmentation axis. Parity 
invariance requires that the probability of forming a given helicity state cannot depend on 
the sign of the helicity j^. However, the fragmentation need not be entirely isotropic; states 
with different magnitudes \j^\ can arise with different probabilities. For the examples we 
will consider here, we can characterize the situation as follows: for a system of light degrees 
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of freedom of spin j, let wj be the probability that fragmentation leads to a state with the 
maximum value of j^. Then the parameter wj may take values between and 1. 

Of the systems available for study, the most interesting is the J = § excited meson 
doublet in the charm system. These states have been identified as the spin-1 Di(2420) 
and the spin-2 (2460), and their decays to the ground state mesons D and D* via pion 
emission have been observed.^) The production of these resonances via fragmentation is 
described by the parameter W3/2, which is the probability of producing the light degrees 
of freedom in a state of helicity ±| along the fragmentation axis, rather than Since 
the Di and D2 peaks are well separated, these two components of the doublet become 
incoherent before they decay, and it becomes appropriate to describe them as a mixed 
state containing Di or D2 with fixed probabilities. If we make the further assumption 
that the initial charm quark is polarized left-handed (it is straightforward to extend this 
to a general or mixed polarization), then the populations of the helicity states of the excited 
mesons are easy to calculate. For example, the helicity states of the D2 are populated with 
the following relative probabilities, where the helicity runs from —2 to +2 across the table: 

p{D*,h)^(^^Ws/2 UI-WS/2) l{l-Ws/2) >3/2 0^. (2.1) 

Given these probabilities, we may now compute the angular distribution for the ob- 
served strong decay D2 ^ D + n. We define 6 to be the angle between the momentum of 
the pion and the fragmentation axis, in the rest frame of the D2- Normalizing to the full 
width r, we find 

1 dV 



Tdcose 



{D; ^ Dtt) = i [1 + 3 cos2 e - 6«;3/2(cos2 e-^)]. (2.2) 



This distribution is independent of the initial polarization state of the charm quark (this 
is not necessarily the case for the decay of the Di; see Ref. 2.). Note that for the isotropic 
production of the light degrees of freedom, W3/2 = |, the pion is also emitted isotropically. 
However, an ARGUS analysis of this decay"^^ found no significant population of the helicity 
states ±2 of the implying that 103/2 is small. Our own fit yields 

W3/2 < 0.24 , 90% conf . (2.3) 

This is an intriguing result, for which we have as yet no satisfactory physical interpretation. 

A similar quantity may soon be measured for the j = 1 heavy charmed baryon doublet, 
the spin- 1 Sc and the spin-| S*. The Ec has been observed with a mass of 2453 MeV; 
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while the S* has not yet been found, its mass is known to be greater than 2530 MeV. In 
any case, the peaks for the two states are known to be well separated. Hence an analysis 
of the fragmentation parameter wi relevant to this system may be performed by observing 
the strong decay E* — > Ag + tt. Defining 9 as before, we find the angular distribution of 
the pion to be 

1 dV 



r d cos 9 



(S* ^ AcTt) = i [1 + 3 cos^ 9 - f tui (cos2 9-1)]. (2.4) 



In this case, isotropic production of the light degrees of freedom corresponds to lUi = |. 

This quantity, when it is obtained, will be useful for another reason. The ground 
state Aq baryon has light degrees of freedom with total angular momentum j = 0; the 
spin of this baryon is carried entirely by the heavy quark. Hence if a sample of Ag 
baryons arises from an initial sample of polarized heavy quarks, they should themselves 
be strongly polarized.^) (We have in mind charm and bottom quarks produced at the Zq 
peak, which have polarizations Pc = 0.67 and Pi, = 0.94, respectively.) The Ag would, 
in fact, carry all of the polarization of the heavy quark if it were produced only directly, 
but there is a potentially significant pollution from Ag's produced via the cascade decays 
(Eg, Eg) Ag + TV. If (as is the case in the charm system) the excited doublet (Eg, Eg) 
lives long enough that spin exchange interactions between the heavy quark and the light 
degrees of freedom have time to take place, the initial polarization of the heavy quark will 
be diluted. The actual polarization of the Ag baryons may be computed in terms of the 
anisotropy parameter wi and a parameter A which describes the probability of producing 
a J = 1 (and isospin one) diquark relative to one with j = (and isospin zero), summed 
over the nine possible helicity and isospin states. If P is the initial polarization of the 
heavy quark, we find the polarization of the sample of Ag's to be 

P^ = i±(i±i^P. (2.5) 

The parameter A is related to one which appears in the Lund Monte Carlo^-* and may be 
estimated to be 0.45, with large errors. If, motivated by the results for charmed mesons, 
we take wi = 0, we find a polarization Pa = 0.72P. Note that the polarization retention 
Pa grows with wi . Finally, the brief treatment given here depends on the assumption that 
the Eg and Eg resonances are well separated. While this is known to be true in the case 
of charm, it may not necessarily be true for bottom. For a more general analysis in which 
this latter possibility is taken into account, see Ref. 2. 
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2.2. Heavy-Heavy Systems 

We now turn to the much rarer possibiHty that the fragmentation process generates an 
additional pair of heavy quarks, rather than just Ught degrees of freedom. Braaten et al.'^^ 
have shown that the production of heavy "onium" states at high energies is dominated by 
such fragmentation processes, which may be calculated perturbatively. A typical process 
would be the production of charmonium from the fragmentation of a high energy charm 
quark produced in Zq decay, cc ^ if) + cc + X . Braaten et al. have shown that such 

a process factorizes to leading order in rric/E., where E is the energy of the il^: 

dViZ'^ ^iP{E) + X) = J2 f dzdt{Zo^i{E/z) + X,iJ,)Di^^{z,iJ,). (2.6) 

i 

Here is the cross-section for the Zq to produce the parton i with energy E/z., and 
Di^tjj (z, jj) is a process-independent fragmentation function for the parton i to fragment 
to a with momentum fraction z. The scale is a factorization scale which is introduced to 
separate the short- and long-distance physics; note that the energy of the V' does not appear 
in Di^^{z, /u), which may be evaluated at a low-energy scale /i ~ rric. The renormalization 
group may then be used to resum the leading logarithms between mc and the high-energy 
scale E. 

The fragmentation functions may be evaluated straightforwardly and perturbatively in 
as(?Ti^), but the results are somewhat messy.^"^-* Of more interest are the fragmentation 
probabilites, obtained by integrating over z. In many cases it can be shown that the 
probabilites P — dz D{z, /j,) are in fact independent of //, at least to leading order. 
Some typical probabilities which Braaten et al. report are* 

P,^^ ^ 1.2 X 10-^ , Pc^^^ « 1.2 X 10-^ , 

Pb^B, ~ 3.8 X 10"^ , Pb^B; ~ 5.4 x 10"^ . ^'^'^'^ 

These results may be extended in a number of ways. First, one may compute not only 
the probability of a charm quark to fragment to a t/^, but the fraction of the time that 
the ip is transversely rather than longitudinally aligned. ®) We define ( to be the ratio of 
the probability of producing a transversely aligned i/j to the total production probability. 
Then a straightforward perturbative calculation gives ( = 0.69, corresponding to a small 
excess of transversely aligned ip^s. (This fraction is independent of the heavy quark mass, 

* I am grateful to T. C. Yuan for providing me with updated versions of these numbers. 
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and is hence the same for T's.) In leading logarithmic approximation, the corresponding 
ratio for gluon fragmentation to V^'s has the same value and is also //-independent. Hence, 
at a hadron collider, where t/^'s are produced by both quark and gluon fragmentation, the 
fraction of transversely aligned ip^s is also C = 0.69. 

This fact has an immediate application to the study of direct '0 production at hadron 
colliders, where it is important to distinguish these ■^'s from those produced via the weak 
decays of b quarks. The alignment of the t/j may be observed in the angular distribution 
of its leptonic decay, i/j — > £~^£~ , which is parameterized as 

, oc 1 + acos^g. (2.8) 
a cos 9 

Here 6 is the angle between either of the leptons and the alignment axis in the ip rest 
frame. For ( = 0.69, we find a = 0.053, a small (5%) asymmetry. By contrast, for the 
weak decay b il>si'^£~ we find a ~ —0.46. Hence the lepton angular distribution may 
provide a useful tool for separating these two sources of t/j^s in a collider environment. 

Finally, we have applied these perturbative methods to estimating the probability 
of fragmentation to baryons which contain two heavy quarks, such as bbq, ccq and bcq 
states.^-* There are a host of possible "doubly heavy" spin-i and spin-| baryons; we found 
the largest production probabilities to be 

P(c^E,„E:J~2xlO-^ 

P(6^Sbe,S*J~5xlO-^ (2.9) 
P(6 ^ Abe) ~ 4 X 10-^ . 

While there is considerable uncertainty in these estimates (they should be trusted at best 
within a factor of two), this is enough to tantalize us with the hope that these states may 
one day be observed at hadron colliders such as an upgraded Tevatron. 

3. Inclusive Heavy Hadron Decays 

There has been much recent progress in the analysis of inclusive decays of heavy 
hadrons. More precisely, these are decays to final states in which some, but not all, of the 
quantum numbers of the decay products are known. For example, we might specify the 
identity and momenta of certain weakly interacting particles, such as leptons or photons, 
or we might restrict the flavor of the final hadrons. The decays are inclusive in that we 
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sum over all final states which can be produced by the long-distance, nonperturbative 
strong interactions, subject to certain constraints which are determined by short-distance, 
perturbative physics. 

We are typically interested in studying quark-level transitions induced by short- 
distance weak interactions, such as 6 — > clv^ b — > S7, or 6 — > se+e". In a lagrangian renor- 
malized far below the weak scale, these interactions appear as local non-renormalizable 
operators with coefficients which contain all dependence on the interesting short-distance 
physics. For example, the low-energy lagrangian contains terms 

^v;,c7p(f - Y')bir{i - 1> (3.f) 

and 

^VtbVZj^mbCrii,) SRa^^hL F^^ , (3.2) 

and we would like to extract the parameters Vd, and C'j{mi,) from experiment. In order 
to do this, however, it is necessary to understand how these quark-level operators induce 
transitions in the physical hadrons which are actually observed. 

The difficulty is that exclusive decays, such as 5 — > K*7 in the case of an underlying 
6 — > 57 transition, are typically not accessible theoretically. This is because they require 
the calculation of nonperturbative hadronic matrix elements such as {K* \ s^a^i^hL for 
which there are no rigorous techniques except in certain situations of enhanced symmetry. 
It is simpler, in a sense, to consider inclusive decays, in which all possble final states are 
summed over. One may constrain the kinematics by fixing the momenta of the nonhadronic 
quanta in the final state, and one may assume in the sum that the fiavor of the final hadron 
may be tagged experimentally. In that case, what one needs instead of sjiCni/bi, \B) 
is the sum over all strange hadrons of the probabilities of producing them, 

W = bLCJ^pSR \X,){Xs\ SRa^^bL \B) . (3.3) 

There is an analogous expression, of course, for inclusive semileptonic weak b decays, and 
so forth. 

If the energy which is released in the decay is large, then the inclusive rate may be 
modeled simply by the decay of a free b quark. The heuristic argument for this is simply 
that the quark-level transition takes place over a time scale which is much shorter than 
the time it takes for the quarks in the final state to materialize as physical hadrons. Hence 
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these two parts of the hadronic decay process do not interfere with each other. Once the 
short-distance interaction has taken place, the probabihty is unity that the quarks which 
have been produced wiU hadronize somehow; we do not need to know the (incalculable) 
probabilities of hadronization to individual final states if we sum over all of them. This free 
quark decay model (FQDM) has been the basis, until recently, of all theoretical studies of 
inclusive heavy quark decays. 

It is intuitively clear that the FQDM will be a good approximation when the energy 
release in the decay is much larger than Aqcd: and when the final quark is far from its 
mass shell after the short-distance interactions have taken place. It has been shown by 
Chay et al}^^ that the FQDM may in fact be justified more rigorously within QCD, via 
an operator product expansion. To do this, one notes that the optical theorem may be 
used to rewrite an inclusive sum such as ( ^.3| ) as the imaginary part of a forward scattering 
amplitude, 

W = 2ImT = 2Im(S|T{6i,cr"^Si?,Si?a^^6L} \B) . (3.4) 

There is now no reference in the expression to any explicit strange hadronic state Xg. Let 
us denote by the momentum of the virtual strange quark in the time-ordered product. 
Since the scale of is set by the available energy m;,, over almost all of the Dalitz plot we 
have P^ ^ ml and the strange quark is far from its mass shell. (The point in the Dalitz 
plot is fixed by the kinematics of the rest of the event, in this case the energy of the photon. 
For semileptonic decays, it would be set by the lepton kinematic invariants.) If we stay 
away from those regions where P^ ~ m^, it is appropriate to expand the nonlocal hadronic 
object T in an operator product expansion, as a series of local operators suppressed by the 
large off-shell momentum. 

The insight of Chay et al. was to notice that it is more convenient to expand directly 
in inverse powers of the large mass nih- This is almost the same thing as expanding in 
the off-shellness of the strange quark, since P^ scales with mf,; it only differs near those 
extreme regions of the Dalitz plot where the expansion breaks down. But this technique is 
considerably more powerful, because it allow us to expand in operators of the heavy quark 
effective theory (HQET). The expansion of T then has the structure 



1 

rrib 



(3.5) 



where the operator is of dimension 2> -\- n and in HQET takes the form 

Or, = cJiTD^^...D^^h. (3.6) 
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The effective fields h{x) which have replaced the quark fields b{x) are obtained by the 
usual HQET rescaling.^) 

The advantage of the HQET formulation is that the heavy quark symmetries may 
be used to compute the necessary matrix elements, up to the desired order in l/rrib- At 
leading order in the expansion of T, the matrix element which is needed is of the form 



This may be computed exactly in HQET; it is proportional to the Isgur-Wise function 
evaluated at the zero-recoil point, where it is absolutely normalized. Chay et al. pointed 
out that the result obtained by truncating T at this leading operator reproduced the free 
quark decay model prediction for the inclusive rate. Hence the FQDM was shown to be the 
first term in a controlled expansion in m^, with all corrections formally of order AqcD/TOfe- 
This was a considerable improvement from its earlier status as an ad hoc, if intuitively 
reasonable, model. 

In fact, Chay et al. showed that the situation is considerably better than it appears 
at first glance, because all corrections of relative order vanish if the decay rate 

is written in terms of the bottom quark, rather than meson, mass. (The distinction is 
significant, since the decay rate is proportional to m^.) The result follows simply from the 
fact that all matrix elements of HQET operators with one derivative vanish at the point 
of zero momentum transfer. This is a surprising and intriguing result, since there are 
certainly l/rrih corrections to the individual exclusive modes; however in the inclusive sum 
they conspire to cancel. 

The recent progress in the field has arisen due to the realization that the leading 
non- vanishing corrections, of relative order may be computed simply in terms of 

two nonperturbative parameters. All the necessary matrix elements of operators O2 may 
be expressed in terms of the matrix elements of two dimension-five operators^^^: 



The parameters Ai and A2 appear also in the expansion of the physical hadron masses in 
terms of the bare quark mass. 



{B\hVh\B) . 



(3.7) 



{B\h{iDfh\B) = 2mBXi, 

{B\ ^-ia^'-)G^,h \B) = QMbM ■ 



(3.8) 



Mb = mb + A- 



1 



(Al + 3A2) , 



Mb* = mt + A- 



2mb 
1 



(Al - A2) . 



(3.9) 



2mb 
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Here A is the energy of the hght degrees of freedom. The physical splitting between the 
vector and pseudoscalar states may be used to extract A2, which represents the additional 
energy of the b quark in the hadron due to its chromomagnetic interactions with the light 
degrees of freedom. We find A2 = 0.12 GeV^. Unfortunately, Ai, which represents the 
b quark's residual kinetic energy in the bound state, is not easily related to any such 
observable. It is probably reasonable to allow it to vary within the range —1 GeV^ < Ai < 
IGeV^ 

We can now obtain the leading nonperturbative corrections to the inclusive width F, 
as well as to other inclusive observables. The expansion will generally take the form 



r oc 



4m^ 



(3.10) 



Here the ellipses denote corrections arising from yet higher order terms in the operator 
product expansion, plus QCD radiative corrections which we have not included. We stress 
that while such radiative effects, such as real and virtual gluon emission, can be as large 
as the terms which we consider here, they may simply be added in as necessary. 

We close by tabulating some recent results which have been obtained with these meth- 
ods. The most extensive work has concerned inclusive semileptonic b ^ c transitions.^^-* 
One may obtain the leading corrections to the full inclusive width, as well as corrections 
to the lepton energy spectrum. Since Prof. Bigi has discussed his own results at this con- 
ference, I have little to add. The corrections to the lepton spectrum are typically a few 
percent, except near the endpoint of maximum energy. There the effect grows dramatically, 
but these divergences are a signal of the breakdown of the operator product expansion and 
are not to be trusted. The computed spectrum must be smeared over lepton energies of 
a few hundreds of MeV for the results to be sensible. Unfortunately it is not possible to 
probe the very endpoint of the lepton energy spectrum in detail using these techniques. 
Hence the study of the charmless semileptonic transition b — > u£u is not practical. 

These techniques have also been applied to the rare 6 — > s transitions B ^ Xs'j and 
B — > Xg£'^£~. For the two-body final state, we find corrections to the inclusive width, ^^'^^^ 



r oc 



and to the average photon energy, -^^^ 



(3.11) 



rub 
2 



1 - 



2ml 



(Ai + 3A2) 



(3.12) 
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Note that it is the quark mass m& which appears in these expressions, not the meson mass. 
The corrections are typically a few percent, which tells us that the FQDM is in fact a 
good approximation when applied to these rare b decays. Finally, we observe that the 
form of the correction to {E^) is just what we would expect, since, according to eq. (|3.9|), 
— (Ai + 3A2)/2mb is precisely the shift in the energy of the free b quark when it is bound 
in the B meson. 

For the Xsi~^i~ final state, the form of the result is too unwieldy to reproduce here, 
and we refer the reader to Ref. 15 for details. We have calculated the leading corrections 
to the lepton invariant mass spectrum, dV/ds, where s = {P£+ + P^-)"^ /ml. They are 
somewhat larger than for b 57, typically on the order of ten percent. However, one 
must be careful to stay away from the tj) resonance region in s, where additional four- 
quark operators can contribute and the operator product expansion is not to be trusted. 
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